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We  survey  the  asymptotic  behavior  of  solutions  to  hyperbolic 
conservation  laws  and  balance  laws,  as  well  as  systems  of  such 
laws  with  the  property  that  each  characteristic  field  is  either 
genuinely  nonlinear  or  linearly  degenerate. 


1.  INTRODUCTION 

We  discuss  the  asymptotic  behavior,  as  t  •+■  of  solutions 
to  the  problem  for  systems  of  strictly  hyperbolic 

balance  laws 


ut  +  f^x  +  8(u,x,t)  =  0 


(1.1) 


in  one  space  variable.  The  vector  field  u(x,t)  takes  values 

in  Rn  and  the  functions  f:  Rn  -*■  Rn  and  g:  Rn  *  R  *  R  Rn 

are  assumed  smooth.  Strict  hyperbolicity  means  that  for  each 

u  £  Rn  the  n  *  n  matrix  Vf(u)  has  n  real  distinct 

eigenvalues  (characteristic  speeds)  A. (u) (u) . 

l  n 

When  the  supply  term  g  vanishes,  (1.1)  reduces  to  the 
system  of  hyperbolic  conservation  laws 


u.  +  f(u) 
t  x 


(1.2) 


Solutions  of  (1.2)  vanishing  at  x  •  ±*  conserve  the  quantity 


u(x,t)dx. 

*  -*00 
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Even  so,  Lp  norms  of  u(x,t),  p  >  1,  may  decay,  as  t 
provided  that  the  solution  is  "spreading  out".  When  f(u)  is 

linear,  characteristic  speeds  are  constant  so  that  each  mode 
travels  with  its  fixed  speed  and  does  not  decay.  When  f(u)  is 
nonlinear,  the  dependence  of  characteristic  speeds  upon  u 
causes  the  convergence  of  characteristics  so  that,  even  if  the 
initial  data  u(x,0)  are  very  smooth,  smooth  solutions  generally 
break  down  in  a  finite  time  and  shock  waves  develop.  Shocks 
and  simple  waves  interact  and  cancel  each  other  out  thus  in¬ 
ducing  decay.  Another  indication  of  decay  is  provided  by 
"entropy  inequalities"  [1] 

H(u)t  +  q(u>x  5  °*  (1*3) 

satisfied  by  admissible  weak  solutions  of  certain  systems  (1.2), 
as  abstract  forms  of  the  second  law  of  thermodynamics.  The 
entropy  decay  mechanism  is  particularly  effective  in  the 
presence  of  strong  shocks. 

Since  in  the  linear  case  there  is  no  dissipation,  one 
should  expect  that  solutions  of  (1.2)  will  decay  faster  when 
f(u)  is  "very  nonlinear".  It  turns  out  that  the  relevant 
condition  on  f(u)  is  genuine  nonlinearity  [2], 

^k^U)  "rk^  ^  u  G  Rn,  (1.4) 

where  r^(u)  is  the  right  eigenvector  of  Vf(u)  associated 

with  A^(u).  Condition  (1.4)  is  not  always  satisfied  by  all 

characteristic  fields  in  systems  of  conservation  laws  arising 
in  continuum  physics.  As  a  matter  of  fact,  there  are  classical 
examples  (e.g.  the  system  of  conservation  laws  of  gas  dynamics) 
where  certain  characteristic  fields  are  linearly  degenerate: 

j  (u)  •  r j  (u)  =0,  u  £  Rn.  (1.5) 

One  should  not  expect  decay  of  components  of  solutions  corre¬ 
sponding  to  such  fields. 

Returning  to  general  balance  laws  (1.1)  we  remark  that 
whenever  the  supply  term  g  is  dissipative  it  collaborates  with 
the  dissipative  mechanism  of  (1.2)  to  speed  up  decay.  On  the 
other  hand,  when  g  is  not  dissipative,  there  is  competition 
with  the  dissipative  mechanism  of  (1.2)  the  outcome  of  which 
depends  on  relative  strength. 


In  the  case  of  a  single  genuinely  nonlinear  conservation  law 
very  precise  information  is  available  which  is  described  in 
Section  2.  In  Sections  3  and  4  we  discuss  asymptotic  behavior 


for  the  single  genuinely  nonlinear  balance  law  and  the.  general 
single  conservation  law.  Known  results  for  systems  are 
surveyed  in  Sections  5  and  6. 


2.  SINGLE  CONSERVATION  LAW.  THE  CONVEX  CASE 

In  this  section  we  discuss  the  asymptotic  behavior  of  solu¬ 
tions  to  the  single  conservation  law 

ut  +  f(u)^  *  f(u)  strictly  convex.  (2.1) 

/ 

In  the  genuinely  nonlinear  case,  f"(u)  >  0,  Lax  [2]  dis¬ 
covered  an  explicit  solution  for  the  initial  value  problem  of 
(2.1)  which  yields  precise  information  on  the  asymptotic  be¬ 
havior  of  solutions.  A  different  approach,  which  also  applies 
when  f"(u)  may  vanish  at  isolated  points,  employs  the  concept 
of  generalized  characteristics  [3]. 

The  following  proposition  [3]  exhibits  the  influence  on  decay 
of  (a)  the  "strength  of  convexity"  of  f(u)  and  (b)  the 
deployment  of  initial  data. 


2.1.  Theorem.  Assume  that 

c|u|P  <  f"(0)  <  C|u|P,  p  >  0,  0  <  c  <  C,  (2.2) 

for  u  in  some  neighborhood  of  0.  Let  u(x,t)  be  an 
admissible  solution  of  (2.1)  on  (-00,00)  *  [ 0 ,<*> )  with 


rx+L 
'  x 


u(y,0)dy  =  0(LS), 


as  L  -*■  °°, 


(2.3) 


uniformly  in  x  £  (-00,00) ,  for  some  s  £  [0,1).  Then 

s-1 

u(x,t)  “  0(tp(1-8)+2  8),  as  t  -*■  “,  (2.4)  . 

uniformly  in  x  £  (-*®,00) . 

2.2.  Corollary.  Assume  that  f(u)  satisfies  (2.2)  in  some 
neighborhood  of  0  and  let  u(x,t)  be  an  admissible  solution 
•of  (2.1)  on  (-“,")  x  [0,«®)  with  u(x,0)  £  L^(-",~), 

*1  <  q  <  oo.  Then  ^ 

u(x,t)  *  0(t  as  t  -*•  (2.5) 


uniformly  in  x  £  (-•,«*). 


The  variation,  of  a  certain  function  of  the  solution  also 
decays,  namely, 

2.3.  Theorem.  Let  u(x,t)  be  an  admissible  solution  of  (2.1) 
on  C-00  ,  00)  x  [O,00)  with  initial  data  u(x,0)  £  L^(-«,").  Then, 
for  any  t  £  (O,00), 

Var{f (u(x,t))  -  u(x,t)f *Cu(x,t))}  <  7  [.  |u(x,0)|dx.  (2.6) 

x  ~  08 

In  connection  to  the  above  theorem,  see  the  discussion 
following  Theorem  4.3,  in  Section  4  below.  ' 


When  the  initial  data  are  periodic,  one  has  information  on 
the  asymptotic  shape  of  the  solution:  * 

2.4.  Theorem.  Let  u(x,t)  be  an  admissible  solution  of  (2.1) 
on  x  tO,00)  with  initial  data  u(x,0)  periodic  of 

mean  u.  Assume  f"(u)  f  0.  Then,  as  t  °°,  u(x,t)  is 
asymptotic  to  order  o(t“^)  to  a  periodic  sawtooth  function. 
The  number  of  "teeth"  per  period  equals  the  number  of  points  in 
a  period  interval  in  which  the  function 


f[u(y, 

J0 


0)  -  n]dx 


attains  its  minimum. 

Interesting  asymptotic  shapes  also  emerge  when  initial  data 
have  compact  support.  In  that  case  solutions  behave 
asymptotically  as  N- waves. 

2.5.  Theorem.  Assume  f"(0)  >  0.  Let  u(x,t)  be  an  ad¬ 
missible  solution  of  (2.1)  on  (-*»,»)  x  [O,00)  with  initial  data 
of  compact  support.  Then 

(0  X  ^  n-(t) 

u(x,t)  -  !  [7  -  f(0)]  +  0(7)  n_(t)  <  X  <  n+(t)  (2.7) 


n+(t)  <  x 


where 


ryt)  -  tf’(0)  +  [2I-tf"(0)]1/2  +  0(1) 


(2.8) 


with  I_, I+  constants  (invariants  of  the  solution)  depending 
upon  the  Initial  data. 


When  f"(0)  vanishes,  the  N-wave  is  deformed,  as  the 
following  example  indicates. 

2.6.  Theorem.  Let  u(x,t)  be  an  admissible  solution  on 
(-00.00)  x  [O,00)  of  the  equation 

ut  +  (u2k)x  =0,  k  =  1,2,...,  1 

with  initial  data  of  compact  support.  Then 


(2.9) 


x  <  n_(t) 


u(x,t)  =4  (^f)2k-1  +  o(t  k)  n  (t)  <  X  <  njt) 


-  v 


h+(t)  <  X 


(2.10) 


where 


2k-l  1 

n-(t)  =  +2k(~Y  I-)  2k  t2k  +  0(1) 


(2.11) 


with  I  ,I+  the  same  invariants  appearing  in  (2.8). 

There  are  also  results  [3,4]  on  the  asymptotic  behavior  of 
solutions  of  (2.1)  when  the  initial  data  approach  different 
constants  u  and  u+  as  x  -+•  and  x  It  turns  out 

that  the  solution  approaches  the  wave  fan  that  solves  the 
Riemann  Problem  [2]  corresponding  to  (u  ,u  ). 

T» 

3.  SINGLE  BALANCE  LAW.  THE  CONVEX  CASE 
We  now  consider  the  single  balance  law 

+  f(u)^  +  g(u,x,t)  *  0*  f(u)  strictly  convex.  (3.1) 

Although  no  analog  of  Lax's  explicit  solution  of  (2.1)  is 
known  for  (3.1),  the  method  of  generalized  characteristics  [5] 
applies  and  yields  precise  information. 

*  As  a  first  example  we  consider  the  case  where  g  is 
flissipative  and  we  generalize  Corollary  2.2. 


3.1.  Theorem.  Assume  that  f(u)  satisfies  (2.2)  in  some 
neighborhood  of  0  and  that 


(3.2) 


ug(u,x,t)  >0,  u  €  x  €  (-®,®),  t  €  [0,®). 

Let  u(x,t)  be  an  admissible  solution  of  (3.1)  on 
(-®,®)  *  [0,®)  with  u(x,0)  £  L*l(-®,®),  1  <  q  <  ®.  Then 

_  1 
p+q+1 

u(x,t)  =*  0(t  ),  as  t  -*• 

uniformly  in  x  £  (-®,00). 


(3.3) 


In  the  next  example  g  is  not  necessarily  dissipative  but 
it  decays  at  infinity  so  that  the  dissipative  mechanism  of  the 
conservation  law  prevails.  < 

3.2.  Theorem.  Assume  that 

|g(u,x,t)|  <  a(x)b(t) ,  u£(-»,®),  x  £  (-®,®),  t  €  ft),”),  (3.4) 

where  a(x)  is  a  bounded  function  such  that  a(x)  -*•  0  as 
|x|  ■*  ®  while  b(t)  €  L^(0,®).  Let  u(x,t)  be  an  admissible 
solution  of  (3.1)  on  (-®,®)  £  [0,®)  with  u(x,0)  6  L*qc  (-“,“) 

such  that  u(x,0)  0  as  |x|  -►  ®.  Then 


u(x,t)  ■  o(l),  as  t  •+■  ®, 
uniformly  in  x  €  (-®,“). 


(3.5) 


It  is  also  possible  to  establish  decay  of  the  variation  of 
solutions.  Here  is  a  typical  result  for  the  periodic  case. 


3.3.  Theorem.  Consider  the  balance  law 


ufc  4-  f(u)x  +  g(u)  *  0 


(3.6) 


with  f"(u)  >  a  >  0,  g(0)  =  0,  and  g'(u)  >  0.  Let  u(x,t)  be 
an  admissible  solution  of  (3.6)  on  (-*,*) ”x  [0,®)  with 
initial  data  periodic  of  period  T.  Then 


Var  u(x,t)  <  , 

[o,t]  -  at 


t  £  (0,®) . 


(3.7) 


•  When  the  initial  data  have  compact  support,  the  solution 
attains  an  asymptotic  profile  analogous  to  the  N-wave  of  the 
conservation  law  (Theorem  2.5).  We  present  the  result  in  the 
context  of  a  simple  concrete  example. 


t 


3.4.  Theorem.  Let  u(x,t)  be  an  admissible  solution  on 


( — 00 , 00 )  x  [0,«)  of  the  balance  law 
2 

Ut  +  (TK  +  u  =  0  (3*8) 

with  initial  data  of  compact  support.  Then 


(  0 


u(x,t)  = 


rz  *  °<i> 


x  +2t 
^  0 


x  <  n  (t) 


n_(t)  <  x  <  n+(t) 


n+(t)  <  x. 


(3.9) 


The  fronts  h  (t),n+(t)  of  the  wave  are  determined 

asymptotically  as  solutions  of  an  equation  which  is  too 
complicated  to  be  illuminating. 


4.  SINGLE  CONSERVATION  LAW.  THE  GENERAL  CASE 


Next  we  consider  the  single  conservation  law 

+  f(u)x  *  0  (4.1) 

without  any  convexity  restrictions  on  f(u).  No  explicit 
representation  of  solutions  is  known.  The  method  of  gener¬ 
alized  characteristics  works  but  not  as  effectively  as  in  the 
convex  case. 


In  general,  decay  is  to  be  expected  when  f"(u)  vanishes 
only  at  isolated  points.  The  rate  of  decay  will  depend  on  the 
"flattness"  of  f(u)  at  points  of  inflexion.  The  following 
result,  established  [6,7]  by  methods  of  topological  dynamics, 
starts  from  very  weak  assumptions  on  f(u)  but,  in  return, 
yields  no  information  whatsoever  on  the  rate  of  decay. 


4.1.  Theorem.  Assume  that  the  set  of  points  on  which  f"(u) 
vanishes  has  no  (finite)  accumulation  point  on  the  real  axis. 
Let  u(x,t)  be  an  admissible  solution  of  (4.1)  on  C-®8,®)  x 
,[0,°°)  with  initial  data  that  are  L  -almost  periodic  on  (-®0,00) 
of  mean  0.  Then  u(x,t)  decays  to  0,  as  t  ■*  «,  in 


In  the  test  case  where  f(u)  has  only  one  inflexion  point, 
one  may  derive  the  asymptotic  behavior  of  solutions  from  the 
convex  case  (Co roller^  2.2)  by  employing  the  ordering  property  of 


admissible  solutions  of  (4.1),  i.e.,  that  u(x,0)  <  u(x,0), 
x  £  (-00,00),  implies  u(x,t)  <  Q(x,t),  x  £  (-00,00),  t  6  [0,") . 

4.2.  Theorem.  Assume  that  f"(u)  >  0  for  u  >  0, 
f"(u)  <  0  for  u  <  0  and 


c|u|P  <  |f"Cu)|  <  C|u|P,  p  >  0,  0  <  c  <  C,  (4.2) 


for  u  in  some  neighborhood  of  0.  Let  u(x,t)  be  an 
admissible  solution  of  (4.1)  on  (-®0,®)  x  [0,*®)  with, 
u(x,0)  £  L<1(-00,00),  1  <  q  <  ®.  Then 

_  1 

u(x,t)  =  0(t  P+q+1),  as  t  -*■  %  ^  (4.3) 

uniformly  in  x  £  (-“,“) . 

Using  a  clever  scaling  argument,  Benilan  and  Crandall  [8] 
obtain  a  decay  estimate  for  orbits  of  Lipschitz  continuous 
semigroups  generated  by  homogeneous  generators.  Since  -f(u) 

1  x 

generates  a  contraction  semigroup  on  L  (-°°,°°),  the  above 
result  applies  to  (4.1),  when  f(u)  is  homogeneous,  and  gives 

4,3,  Theorem.  Assume  that  for  some  u  >  2 

f(Xu)  =  Aaf(u),  u  £  (-“>,“>),  Xe[0,~).  (4.4) 

Let  u(x,t)  be  an  admissible  solution  of  (4.1)  on  (-®,00)  x 
[0,°°)  with  u  (x ,  0)  €  (-00,00) .  Then 

Var  f(u(x,t))  <  I”  |u(x,0)  |dx,  t  €((),<*).  (4.5) 

In  particular,  (4.5)  implies  |u(x,t)|a  *  0(t  1)  which  is 
in  agreement  with  (2.5),  (4.3).  Note  that  (4.5)  also  yields 
information  on  the  decay  of  u  (x,t).  In  the  convex  case, 

i  ia  C 

f(u)  *  |u|  ,  (4.5)  reduces  to  (2.6).  It  is  not  known  whether 
(2.6),  or  any  analog  of  it,  holds  for  general  conservation  laws. 

s  The  decay  of  solutions  with  periodic  initial  data  was  in- 
vesitgated  by  Greenberg  and  Tong  [9]  who  obtain  the  following  - 

4.4.  Theorem.  Assume  that  f(u)  is  as  in  Theorem  4.2  and  let 
u(x,t)  be  an  admissible  solution  of  (4.1)  on  (-*,00)  x  [0,®) 
with  initial  data  periodic  of  mean  0.  Then  there  is  a  sequence 


/ 


{t  },  with  t  -*■  ».  such  that 
n  n 

_  1 

u(x,tn)  »  0(tn  P+1).  n  °°,  (4.6) 

uniformly  in  x  6 

The  asymptotic  shape  of  periodic  solutions,  under  the 
assumptions  of  Theorem  4.4,  has  also  been  studied  [9,10]. 

Our  understanding  of  the  asymptotic  behavior  of  solutions 
to  the  general  conservation  law  (4.1)  is  still  imperfect. 
Additional  research  will  be  required  for  the  completion  of  the 
program. 

* 


5.  SYSTEMS  OF  CONSERVATION  LAWS 

All  available  information  on  the  asymptotic  behavior  of 
solutions  to  systems  of  conservation  laws  has  been  obtained 
through  explicit  studies  of  wave  interactions,  within  the  frame¬ 
work  of  the  construction  scheme  of  Glimm  [ll],  and  is  thus  re¬ 
stricted  to  solutions  of  small  oscillation  or  variation. 

We  discuss  first  the  case  where  the  system  is  genuinely 
nonlinear,  i.e.,  (1.4)  is  satisfied  for  k  =  l,...,n.  The 
crucial  observation  is  that  weak  waves  propagate  and  interact  as 
waves  of  a  single  conservation  law,  modulo  an  "error"  which  is 
at  most  of  second  order  in  the  strength  of  interacting  waves. 

The  general  strategy  is  to  show  that  the  cumulative  effect  of 
these  errors  is  insignificant  and  thus  decay#  of  waves  in  systems 
is  governed  by  precisely  the  same  laws  that  yule  decay  of  waves 
in  a  single  conservation  law.  Systems  of  two  equations  are 
special  in  that,  in  Riemann  invariants  coordinates,  interaction 
errors  are  actually  of  third  order  and  this  makes  the  study  of 
decay  easier.  The  earliest  result  is  due  to  Glimm  and  Lax  [12]: 

5.1.  Theorem.  When  the  initial  data  u(x,0)  have  small 
oscillation,  there  is  a  solution  u(x,t)  on  (-*9,°°)  *  [0,<»)  to 
the  strictly  hyperbolic,  genuninely  nonlinear  system  (1.2)  of 
two  equations.  Furthermore,  if  u(x,0)  is  periodic  of  mean  0, 

•  u(x,t)  -  0(7),  as  t  -*>  »,  (5.1) 

while  if  u(x,0)  has  compact  support, 

-  I 

Var  u(x,t)  =  0(_t  ^),  as  t  -*•  00 . 

x 


(5.2) 


Starting  from  the  above  result,  DiPerna  [13]  was  able  to 
show  that  in  systems  of  two  equations  the  two  characteristic 
fields  decouple  at  a  rate  0(t“3/2)  and  each  one  approaches, 
at  a  rate  0(t-^-/^),  an  individual  N-wave  of  the  type  exhibited 
in  Theorem  2.5.  After  a  careful  analysis,  Liu  [14]  observed  that 
in  a  genuninely  nonlinear  system  of  n  equations  interactions 
between  waves  of  different  families,  which  produce  errors  of 
second  order,  are  rapidly  completed,  since  these  fields 
propagate  with  distinct  characteristic  speeds.  On  the  other 
hand,  interactions  between  waves  of  the  same  family,  which  per¬ 
sist  longer,  produce  errors  of  the  third  order.  He  thus 
establishes  the  following  , 


5.2.  Theorem.  When  the  initial  data  u(x,0)  have  compact 
support  and  small  total  variation,  there  is  a  solution 
u(x,t)  on  (-00,00)  x  [0,°°)  to  the  strictly  hyperbolic,  genuinely 
nonlinear  system  (1.2)  of  n  equations  and 


1 

~  9 

Var  u(x,t)  =  0(t  ),  as  t  -*■ 

x 


(5.3) 


Furthermore,  as  t  +  ”,  u(x,t)  approaches  in  L  at  a  rate 
0(t-l/6),  a  system  of  n  N-waves  propagating  with  speeds 

MO),. ...A  (0). 

1  n 


It  is  not  known  whether  the  rate  0(t 
theorem  is  optimal. 


-1/6 


)  in  the  above 


Next,  we  consider  systems  with  the  property  that  each 
characteristic  field  is  either  genuninely  nonlinear  or  linearly 
degenerate.  DiPerna  [15]  studied  the  asymptotic  behavior  of 
solutions  to  systems  in  the  above  class  that.,  are  endowed  with  an 
entropy.  Subsequently,  Liu  [16]  established  more  precise  re¬ 
sults  based  exclusively  upon  considerations  on  wave  interactions. 


5.3.  Theorem.  Consider  a  strictly  hyperbolic  system  of  n 
equations  each  characteristic  field  of  which  is  either  genuinely 
nonlinear  or  linearly  degenerate.  When  the  initial  data  have 
compact  support  and  small  total  variation,  a  global  solution  to 
the  Cauchy  problem  exists  and,  as  t  °°,  each  genuinely  non¬ 
linear  characteristic  field  approaches  in  L^-,  at  a  rate 
0(t-l/6),  an  N-wave,  while  every  linearly  degenerate  character¬ 
istic  field  approaches  in  L  ,  at  a  rate  0(t-1'2),  a  travelling 
Vave  (generalized  contact). 

Liu  also  considers  [16,17,18]  initial  data  which  approach 
different  constant  states  u_  and  u+,  as  x  -*■  -00  and  x  -*■  00 , 
and  shows  that  the  resulting  solution  approaches  the  wave  fan 


that  solves  the  Riemann  problem  associated  with  (u_,u+).  Thus, 
the  investigation  of  asymptotic  behavior  within  the  framework  of 
small  solutions  and  Glimm's  method  is  virtually  complete.  In 
contrast,  very  little  is  known  about  the  asymptotic  behavior  of 
solutions  with  large  initial  data  or  of  solutions  to  systems 
that  are  not  genuinely  nonlinear. 


6.  SYSTEMS  OF  BALANCE  LAWS 

The  investigation  of  the  asymptotic  behavior  of  solutions 
to  systems  of  balance  laws  (1.1)  has  only  begun  recently.  Using 
an  adaptation  of  the  method  of  Glimm,  Liu  [19]  studies  systems 
of  balance  laws 

ut  +  f(u)^  +  g(u,x)  =0  *  (6.1) 

under  the  assumption  that  g  decays  rapidly  as  |x|  ■+•  00  and  no 
characteristic  speed  is  0.  Liu's  scheme  is  based  upon  an 
approximate  resolution  of  discontinuities  by  (6.1),  analogous 
to  the  Riemann  problem  solution  for  (1.2)  which  is  the  building 
block  of  Glimm's  scheme.  A  typical  result  is  that  when  the 
initial  data  u(x,0)  have  small  total  variation  and  approach 
constant  states  u  and  u+  as  x  -*■  and  x  -*■  +°°,  there  is 
a  solution  u(x,t)  of  (6.1)  on  (-°°,“)  x  [0,°°)  which 
approaches,  as  t  00 ,  a  fan  of  waves,  which  do  not  produce 
interactions  and  cancellations  and  essentially  consist  of  simple 
and  shock  waves  of  the  associated  conservation  law,  together  with 
a  steady  state  solution 

+  g(u,x)  =  0  (6.2) 

♦ 

of  (6.1)  which  occupies  the  center. 

It  is  not  known  what  happens  in  the  "resonance"  situation 
where  one  of  the  characteristic  speeds  is  0.  Clearly,  a  lot  of 
work  is  still  needed  in  order  to  understand  completely 
asymptotic  behavior  in  systems  of  balance  laws. 
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